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Studies of the ultrarelativistic collisions of hadrons and nuclei at different centrality and energy 
enable to explore the QCD phase diagram in a wide range of temperature and baryon density. 
Long-range correlation studies are considered as a tool, sensitive to the observation of phase 
transition and the critical point. In the present work, a Monte Carlo model of proton-proton, 
proton-nucleus, and nucleus-nucleus collisions is applied to heavy and light ion collisions at the 
cms energy range from a few up to several hundred GeV per nucleon. The model describes the 
nuclear collisions at the partonic level through interaction of color dipoles and takes into account 
the effects of string fusion, which can be considered as an alternative to relativistic hydrodynam¬ 
ics way of describing the collective phenomena in heavy-ion collisions. The implementing of 
both the string fusion and the finite rapidity length of strings allowed to consider the particle pro¬ 
duction at non-zero baryochemical potential. We calculated the long-range correlation functions 
and correlation coefficients between multiplicities and transverse momentum at several energies 
for different colliding systems and obtained predictions for the experiment. 
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1. Introduction 

The exploration of the phase diagram of strongly interacting matter and a search for the onset 
of deconfinement and the critical point is one of the main subjects of heavy ion physics. There 
is a general consensus that at zero baryochemical potential the deconfinement phase transition is 
a smooth cross-over [1], which is mainly based on the lattice QCD calculations. However, the 
numerical lattice study of the behaviour of the QCD matter at non-zero baryon density, due to 
the sign problem, is not straightforward [2, 3, 4]. According to existing calculations [5, 6, 7] 
it is expected that at high net baryon density, this transition is of the type of the first order, which 
suggests an existence of a critical point in the phase diagram at the intermediate net-baryon density. 
However, in several works a smooth crossover has been obtained at any value of baryochemical 
potential [8, 9]. 

The experimental investigations on the QCD phase diagram are related to the study of the col¬ 
lisions of the ions at high energy [10, 1 1]. The study of the QCD phase diagram is a part of physical 
program of the NA61 experiment at the SPS, experiments at RHIC and also future detectors CBM 
at FAIR and MPD at NICA [12, 13, 14, 15]. Among the observables, being used in these studies, 
the most sensitive ones, such as a collective flow, correlations, fluctuations, require event-by-event 
analysis [16, 17]. Particularly, the studies of long-range correlations between variables taken from 
two different rapidity windows are included in the research program of the experiment NA61 as a 
tool, sensitive to the observation of phase transition and the critical point [17]. 

For the correct interpretation of experimental findings and evaluation of the sensitivity of ex¬ 
perimental methods, the theoretical modelling of the evolution of a heavy-ion collision with the 
explicit calculation of observables in the conditions, close to the experimental ones, is required. 
Due to the complexity of the evolution process of a heavy-ion collision and inapplicability of QCD 
perturbation theory in the low momentum region, semiphenomenological models are widely ap¬ 
plied in this area. 

One of the models, which is used to describe the properties of the initial dense state of the 
strongly interacting matter that occurs immediately after the nuclear interaction of high energy, is 
the model for the formation and fusion of quark-gluon strings [18, 19, 20, 21, 22]. According to 
this model, the hadrons, produced in the soft region, can be described as a result of decay of colour 
field tubes - strings, which are formed between the interacting partons. With increasing energy 
and the mass number of colliding nuclei, the density of strings grows, and they begin to overlap 
forming an area in a transverse plane with stronger color field ("clusters"). In the limit of high 
density, it is expected that all cross-section plane will be a single cluster, which is supposed to be in 
its physical properties equivalent to the quark-gluon plasma. Convincing evidence in favour of the 
fusion of the strings are discovered experimentally growth of the mean transverse momentum with 
multiplicity in hadron and nuclear collisions and a significant decrease in the yield of multiplicity 
in the collisions of heavy ions compared to models with independent strings, which is compatible 
with the existing experimental data on multiplicity. It is assumed that fusion and percolation of 
strings responsible for the appearance of the ridge in two-particle correlations [23, 24]. In the 
recent papers [25, 26, 27] it was shown that the equation of state of QGP at zero chemical potential 
obtained in the colour string percolation model is in an excellent agreement with the lattice results. 
In addition, the approach has been successfully applied for the determination of the shear viscosity 
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[28] of QCD matter, with was found in a good agreement with experimental data in a wide range 
of temperatures. 

In the framework of the string fusion approach, the critical behaviour is expected when the 
processes of string fusion and percolation come into play, what can be considered as a possible 
way of quark-gluon plasma formation. Around percolation threshold, strong fluctuations in colors 
of strings appear what lead to large fluctuations in some observables, which one can find by the 
event by event analysis [17]. In the present work, the model is applied to heavy and light ion 
collisions at the center of mass energy range from a few up to several hundred GeV per nucleon. 
In the next sections, we very briefly describe the main features of the Monte Carlo model with 
string fusion and explain how it can be applied at non-zero baryon density. Then, we introduce 
the observables (correlation functions and correlation coefficients) used in the study. Finally, we 
present the results and give our conclusions. 

2. The Monte Carlo model 

The present model [29, 30] describes nuclear collisions on the partonic level, without referring 
to the Glauber picture of the independent collisions of the nucleons. Initially, the nucleons are 
supposed to be distributed in a nucleus according to the charge density distribution p(r). We 
used harmonic oscillator model for 7 Be and 9 Be nuclei and Woods Saxon distribution for heavier 
species [31]. 


Nucleus 

Model 

pW 

Parameters 

2 °8p b 

Woods Saxon 

Po 

1 +exp [{r-R)/d\ 

R = 6.63 fm, d = 0.545 fm 

197 Au 

Woods Saxon 

Po 

1 +exp [(r-R)/d\ 

R = 6.38 fm, d 0.535 fm 

40 Ca, 40 Ar 

Woods Saxon 

Po 

1 +exp [(r-R)/d\ 

R = 3.53 fm, d = 0.542 fm 

7 Be 

Harmonic oscillator 

po (1 + a(r/r 0 ) 2 ) exp (- (r/r 0 ) 2 ) 

r 0 =1.77fm, a=0.327 

9 Be 

Harmonic oscillator 

po (1 + a(r/r 0 ) 2 ) exp (- (r/r 0 ) 2 ) 

r 0 = 1.791 fm, a=0.611 


Table 1: Nuclear distributions and parameters used in the present study. 


In the model, each nucleon contains a valence quark-diquark pair and several sea quark- 
antiquark pairs. The number of sea pairs is generated according to Poisson distribution around 
some mean value, which is energy-dependent and is adjusted to describe the total inelastic cross 
section in pp collisions. The transverse coordinates of these partons are distributed around the cen¬ 
ter of each nucleon according to two-dimensional Gauss distribution with a mean-square radius tq. 
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It is assumed that quark-diquark and quark-antiquark pairs form dipoles. An elementary interaction 
according to the color dipole approach [32, 33] is realized with the probability amplitude given by: 
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where Kq is modified Bessel function of the second order, ( f\ frf), {f\ , r' 2 ) are transverse coordinates 
of the projectile and target dipoles, and a s - effective coupling constant, r max is characteristic 
confinement scale. According to this formula, two dipoles interact with higher probability, if the 
ends of dipoles are close to each other in the transverse plane, and (other equal) if they are wide. 

In our Monte Carlo model it is assumed, that if there is a collision between two dipoles, two 
quark-gluon strings are stretched between the ends of the dipoles, and the process of string frag¬ 
mentation gives observable particles. The particle production of a string is assumed to go uniformly 
between the string rapidity ends y m i n and y max [29, 34], which are related to the longitudinal mo¬ 
mentum fraction, carried by the parons forming the string. The emission of the charged particles is 
generated with mean number per rapidity (p), independently in each rapidity interval, with Poisson 
distribution. 

The transverse position of string is assigned to the arithmetic mean of the transverse coordi¬ 
nates of the partons at the ends of the string. Due to finite transverse size of the strings they overlap, 
that in the framework of string fusion model [18, 19] gives the source of a higher tension. We used 
the lattice variant of the model [35, 36], according to which strings are supposed to be fused if 
their centres occupy the same cell. The area of each cell is equal to the transverse string area 
(S str = Tnftr)- I n calculations, we also assume local string fusion scenario [35, 36]. Each bin in 
rapidity with an integer number of overlapped strings is processed separately. Mean multiplicity of 
charged particles and mean p t originated from the cell where k strings are overlapping are obtained 
according the following: 

(A0* = roVic, (p t ) k = p Q Vk. 


Here po and po are mean charged multiplicity from one single string per rapidity unit and mean 
transverse momentum from one single string. The relations allows to calculate long-range correla¬ 
tion functions and correlation coefficients between multiplicities and mean transverse momentum 
of charged particles [30]. Important, that po and po cancels out in the correlation observables 
defined below and do not influence the final result on correlation coefficients, which makes the 
calculations more robust. 

The parameters of the model a s , ro, r max , po are supposed to be energy-independent and have 
been constrained from charged particles multiplicity in a wide energy range (from ISR to LHC) 
in pp collisions and in minimum bias p-Pb and Pb-Pb collisions at the LHC energy. The detailed 
procedure of parameters tuning is described in [34]. 


Important to note, that the implementation of the string fusion in the model, described above, 
allows us to separate phenomena taking place (in the same event) in different rapidity regions. At 
the SPS energies, the contribution of sea quarks and corresponding strings is not significant. There¬ 
fore, the valence strings dominate. Due to the fact that the diquark carries out larger longitudinal 
momentum fraction [29], the valence string, formed by a quark-diquark pair, is asymmetric. Most 
of the particles in forward rapidity region are coming from valence strings from their diquark parts. 
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The approximate particle content of the valence rapidity string is shown in Fig. 1. The diquark 
end of the string, which is generally located in the forward rapidity region, is characterized by 
large baryon yield. In contrary, the opposite end, which carries out a quark, emits mostly pions, 
which happens around midrapidity. Altogether, it provides a distribution of a net-baryon density 
over rapidity with large baryon redundancy in the forward region, which has been obtained in the 
experimental data (Fig. 2). We note, that the valence string model, developed in [37], successfully 
describes net baryon yields in a wide energy range. 



Figure 1: A sketch of a particle composition, coming from one string: fractions of protons, antiprotons and 
mesons as a function of rapidity. 



Figure 2: Net baryon rapidity distribution in Au+Au collisions. Points - experimental data at RHIC, lines - 
results of calculations in a valence string model [37]. 
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3. Long-range correlations 

The long-range correlations between observables in two separated rapidity windows (back¬ 
ward B and forward F) are numerically characterized by correlation functions and correlation co¬ 
efficients. The correlation function is defined as an average of the observable on backward window 
at fixed value of the observable in the forward window: 


m=( b) f . 


(3.1) 


Both B and F could be either A C h - the number of charged particles in the rapidity window or 
p t - mean in the event transverse momentum of charged particles in the given window: 


] N ch 


(3.2) 


Correspondingly, tree types of correlation functions are defined [35, 36]: n-n (/ n _ n ), 
pt-n (/p t —n) and pt-pt (/ Pt - Pt ). 


fn-n(np) — (' n B] nF i fp t -n(np) — {Pt B )n F ’ /pt-pt ( n ^) — (ftfi) PtF - (3-3) 

The correlation coefficient represents the slope of the correlation function: 

. df{F) 

VS-F = — 7 — F=<F>- (3.4) 

dF 

In order to reduce experimental and theoretical uncertainties, it is often useful to switch to 
normalized variables: B -+ B/(B),F -+ F/(F). In this case, the p t — n correlation coefficient 
becomes dimensionless, and also both n — n and p t — p t correlation coefficients do not change in 
case of symmetrical rapidity windows. In normalized variables: 


< nf > d <hb> 


<n B > 

dn F '" F=< " F> ’ 

< Up > 

d < PtB > 


< PtB > 

dnp 

\np=<np>^ 

< PtF > 

d <p tB > 


<PtB > 

dpt F 

\PtF = <PtF> 


(3.5) 

(3.6) 

(3.7) 


There is an alternative definition on long-range correlation coefficient [38]: b corr = - 

The definitions are equivalent in case of linear correlation functions, and both of them are used in 
theoretical and experimental studies [39, 40, 41, 42]. 

In our calculations we took p+p, 7 Be+ 9 Be, p+Pb, Ar+Ca, Au+Au and Pb+Pb collisions at the 
colliding energies y/s from 5 GeV to 62.4 GeV. The set of the colliding systems and energy range 
have been selected in accordance with the experimental program of NA61 Collaboration, BESII 
program at RHIC and future detectors at FAIR and NICA facilities. 

For p+p, Be+Be, p+Pb collisions minimum bias events have been taken; for heavier species 
(Ar+Ca, Au+Au, Pb+Pb) two centrality classes have been considered: the central events, which 
correspond to the number of participant nucleons Apart > (A + B)/2, and the peripheral ones, with 
Apart < (A + B)/2. Here, A and B are mass numbers of the colliding nuclei. 
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We selected three rapidity windows configurations: 

(-1 ; 0 ) - (0 ; 1 ); 

( 0 ; 1 ) - (1 ; 2 ); 

( 0 ; 1) - (2 ; 3). 

These windows correspond to mid-rapidity configuration (first two) and middle-forward case (third 
one). The forward window in the third configuration probes the rapidity region, where significant 
net baryon density is expected (see Fig. 2). 

4. Results 

4.1 Correlation functions 

Fig. 3 shows all three types of correlation functions. The n-n correlation function is obtained 
to be close to linear at low values of up, which correspond to peripheral collisions. The slope of 
the correlation function decreases with np and tends to saturate at high values of np. The results 
on pt-n correlations demonstrate small positive pt-n slope of pt-n correlation function. The pt- 
pt correlation function is found strongly non-linear in minimum bias collisions. The predictions 
are similar to the ones, obtained for Pb+Pb collisions at the SPS energy in string fusion approach 
[35, 36] and experimental data of the NA49 experiment [43]. 
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Figure 3: Correlation functions in Ar+CA collisions at y/sNN=17GeV: n-n correlation (left), pt-n correlation 
(right) and pt-pt correlation (bottom). Rapidity windows configuration is (-1;0), (0;1). 
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4.2 Correlation coefficients 


In the Fig. 4 the multiplicity-multiplicity correlation coefficient in proton-proton collisions 
is plotted as a function of the center-of-mass collision energy per nucleon yAyy. The results 
demonstrate monotonic increase of the Z? n _ n with y/s. 


bn—n 




— (-1;0), (0;1) 
-®- (0; l), (i;2) 
- (0; 1), (2;3) 



Figure 4: Dependence of n-n correlation coefficient on the collision energy in proton-proton, Be+Be and 
p+Pb collisions. Three rapidity windows configurations are represented by different colors. 


bn—n 




Figure 5: Dependence of n-n correlation coefficient in three rapidity windows configurations on the collision 
energy in Ar+Ca and Au+Au collisions. The color notation is the same as in Fig. 4. 
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In general, the n-n correlation coefficient is higher in mid-rapidity and decrease going from 
mid-rapidity to middle-forward configuration. Such behaviour is in agreement with the previous 
studies [30, 44, 45], where it was shown that Z? n _ n decreases with an increase of the rapidity gap 
and going to the forward rapidity region. 

The results on n-n correlations in proton-ion and light ion collisions are also shown in Fig. 4. 
The correlation coefficient is significantly higher than the one in pp collisions, but no saturation 
with the energy is achieved. 

Multiplicity-multiplicity correlation coefficients for heavier ion species are plotted in Fig. 5. 
The results are split into two different centrality classes, as described above. The Pb+Pb collisions, 
which have a similar pattern as Au+Au, are not shown. It is found, that the b n - n is less in the 
central class than in the peripheral one. Such behaviour is in a correspondence with the string 
fusion predictions [44]. This fact is also in agreement with a saturation trend of the correlation 
function (Fig. 3). However, one should keep in mind that this observable is subject to the centrality 
class width dependence [44], and the central class in our calculations seems to be narrower. So, in 
order to provide the more convincing predictions, the calculations in narrow centrality classes are 
required. 


b P,-n 



sion energy in various colliding systems. The color notation is the same as in Fig. 4. 
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The results on pt-n correlation coefficient are shown in Fig. 6. We find the smooth monotonic 
behaviour of b vt - n with the collision energy, similar to the one, observed for multiplicity corre¬ 
lations. We should note, that at minimal energy, the pt-n correlation is practically zero, which is 
not true for n-n correlations. It may be related to the fact that n-n correlations are present also in 
the absence of any collectivity, due to the fluctuations in the number of the emitting sources [46]. 
In contrary, the pt-n correlations appear only in the presence of the string fusion, and their sign 
depends on the variance of the initial sources [47]. We should note that the magnitude of the pt-n 
correlations increases significantly from light to heavier nuclei. However, this observable also suf¬ 
fering from the dependency on the centrality selection [44], so a more deep centrality investigation 
would improve the calculations. 


bpt-pt 



sion energy in various colliding systems. The color notation is the same as in Fig. 4. 


Fig. 7 shows the energy dependence of the transverse momentum correlation coefficient. Con¬ 
trary to the previously discussed observables, Z? pt - p t is a correlation between two intensive variables 
[36], and it is more robust towards the centrality selection, because it is less influenced by the vol¬ 
ume fluctuations due to variation in participant nucleons number and impact parameter. 
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In the results of the calculation of fo pt _ p t one observes a monotonic increase with the collision 
energy mostly. However, in Ar+Ca there is a non-monotonic dependence with a minimum around 
y/s= 17 GeV. One should stress that the pt-pt correlation coefficient behaves non-monotonically 
only in case of middle-forward rapidity configurations, while in the other two cases (midrapidity 
ones) the dependence is smooth. This picture suggests the existence of a phase transition in the area 
of forward rapidity (note that the larger rapidity corresponds to higher baryon chemical potential) 
with smooth crossover in midrapidity (where the baryon density is minimal). Note that similar non¬ 
monotonic behavior has been obtained experimentally in the energy range 10-20 GeV for different 
fluctuational observables (for example, the deviations from Poisson of net-proton fluctuations at 
RHIC [48]), and it was attributed to the existence of the phase transition and the critical point on 
the QCD phase diagram. In the Fig. 8 we show a more detailed study of pt-pt correlations in Ar+Ca 
collisions in the energy range between 5 and 20 GeV. 


V-pt 




Vs, GeV Vs, GeV 

Figure 8: Dependence of pt-pt correlation coefficient in three rapidity windows configurations on the colli¬ 
sion energy in Ar+Ca collisions. 

We performed these addition calculations in order to study the shape of the observed dip in pt- 
pt correlation energy dependence. The results confirmed previous calculation. The central position 
of the dip in peripheral collisions is found at higher energy than in central, which suggests that the 
phase transition happens earlier in more central collisions. Overall, the obtained picture advises, 
that being applied to the range of non-zero baryon chemical potential, the present model, based 
on the string fusion approach, reproduces the phase transition one expects in the QCD equation of 
state. 

The calculations highlight the importance of Ar+Ca run in the beam energy and system size 
scan by the NA61 Collaboration at the SPS scheduled on 2015 year [49]. The experimental study 
of the long-range pt-pt correlations, as well as different strongly intensive observables [50, 51], 
would be crucial for the improvement of our understanding of the phase diagram of the strongly 
interacting matter. 
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5. Conclusions 

String fusion approach to the quark-gluon plasma formation at non-zero baryon chemical po¬ 
tential has been proposed. A model for the string fusion, accounting finite rapidity width of strings, 
for pp, pA and AA collisions is developed and applied at the SPS energies. Long-range correlation 
functions and coefficients are studied. Smooth monotonic behavior of n-n and pt-n correlation with 
energy with non-monotonic pt-pt correlation dependence in Ar+Ca collisions have been obtained. 

A more detailed scan, including calculation of correlations in narrow centrality classes, is required, 
as well as an extension of the model for net-charge and net-proton fluctuation and correlation stud¬ 
ies, exploring the strongly intensive variables. 
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